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We present a sum-rule extraction of the decay constants of the charmed mesons D and Ds from the 
two-point correlator of pseudoscalar currents. First, we compare the perturbative expansion for the 
correlator and the decay constant performed in terms of the pole and the running MS masses of the 
charm quark. The perturbative expansion in terms of the pole mass shows no signs of convergence 
whereas reorganizing this very expansion in terms of the MS mass leads to a distinct hierarchy of the 
perturbative expansion. Furthermore, the decay constants extracted from the pole-mass correlator 
turn out to be considerably smaller than those obtained by means of the MS-mass correlator. Second, 
making use of the OPE in terms of the MS mass, we determine the decay constants of both D and 
, Ds mesons with an emphasis on the uncertainties in these quantities related both to the input QCD 

parameters and to the limited accuracy of the method of sum rules. 
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PACS numbers: 11.55.Hx, 12.38.Lg, 03.65.Ge 



INTRODUCTION 



^ , The extraction of the decay constants of ground-state heavy pseudoscalar mesons within the method of QCD sum 
' rules P, Q poses a complicated problem. First, one derives an operator product expansion (OPE) for the correlation 
, function 



Hip') = ^ J d'xe^P^O\T [Mx)jI{0)) |0) (1.1) 
of two pseudoscalar heavy-light currents 

^5(2;) = (mg + m)q{x)ij5Q{x). (1.2) 

Second, one considers the sum rule for this correlator. The sum rule is nothing but the expression of the fact that the 
representation of the Borelized correlator (|l.ip . n(p^) 11{t), in the language of the intermediate hadron states is 
equal to the OPE for this correlator: 

00 cx^ 

n(r) - /2 Af^e-A^Q- + J dse-^>i,adr(s) = J ds e-'^ p^.^tis, fi) +U„{t, fi). (1.3) 

Sphys (mQ+niy 

Here, Mq denotes the mass of the pseudoscalar meson Pq containing the heavy quark Q while /g is its decay constant: 

(rng + to)(0|q175Q|Pq) = /qA^§. (1.4) 

For the correlator (|l.ip . Sphys = {Mvq + M^^)^ is the physical continuum threshold, Mvq being the mass of the vector 
meson containing Q. Obviously, for large values of r the contribution of the excited states decreases faster than the 
ground-state contribution and therefore TI{t) is dominated by the ground state. 

The perturbative spectral density is obtained in the form of an expansion in terms of the strong coupling as(/i): 

Ppe.(., A^) = P^'\s) + ^p^^Hs) + (^Ypi^Hs) + ■■■. (1.5) 

TT \ TT / 

Clearly, the correlator (jl.ip does not depend on the renormalization scale fi; however, both the perturbative expansion 
truncated at fixed order in ag and the truncated power corrections npowcr(''": m) given in terms of the condensates and 
the radiative corrections to the latter depend on fi. Moreover, the relative magnitudes of the lowest-order contributions 
strongly depend on the choice of the renormalization scheme/scale. 
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Unfortunately, the truncated OPE allows one to calculate the correlator only at not sufficiently large r, such that 
the excited states give a sizable contribution to n(T) in the corresponding r-range. In principle, the physical spectral 
density above the threshold might be measured experimentally; in practice, however, it is unknown. Therefore, one 
adopts the concept of duality to relate the contribution of the excited hadron states to the perturbative contribution: 
perturbative QCD spectral density Pport(s) and hadron spectral density /3hadr(s) are close to each other at large values 
of s; thus, for sufficiently large values of the parameter s, (far) above the resonance region, one has the duality relation 

oo oo 

y dse"*'^Phadr(s) = J dse~'''^ppcit{s). (1.6) 

In order to express the excited-state contribution by the perturbative contribution, we need to extend this relationship 
down to the value of the hadronic threshold Sphys- However, one has to be careful: the spectral densities Ppcrt('S) and 
/5hadr(s) are obviously different in the region near Sphys- Therefore, one finds 

oo 

J dse"^'^Phadr(s) 

where Seff(T) is different from the physical threshold Sphys- A crucial (albeit rather obvious) observation is that, for the 
same reason which causes Sesir) ^ Sphysi ■Sefi(T) has to be a function of the parameter r to render relation (|1.7p exact. 
By virtue of (|1.7p we may rewrite the sum rule (|1.3p as 

Sott(T) 

fi^M^e-^'i^ = J dse-^>pcrt(s,/x)+npowcr(T,Ai) EEndual(r,Scff(T)). (1.8) 

(mQ+m)2 

We refer to the right-hand side of this equation as the dual correlator. Evidently, even if the QCD inputs ppcrt('S,/^) 
and Ilpowcr (■'■, are well-known, the extraction of the decay constant requires a further criterion for fixing the effective 
continuum threshold Scff (''')• 

Noteworthy, Eq. ()1.8|) offers another way to convince oneself that Soff(r) must be a function of r. In fact, the log slope 
on the left-hand side of (jl.Sp is independent of r and is equal to Mq (which may be exactly known from experiment). 
Consequently, to guarantee the same r-behaviour on the right-hand side of (jl.Sp . the effective threshold must be, in 
general, a function of r. In the literature the approximation of the threshold by some constant sq independent of r is 
widely used. The corresponding dual correlator, IlduaK''', sq), should therefore lead to the presence of a contamination 
of excited states on the left-hand side of Eq. ()1.8p . In principle, one may develop models for excited states in order to 
estimate (and subsequently remove) such a contamination. It is, however, clear that ultimately the same effect can be 
equivalently reached by considering an explicit r-dependence of the effective continuum threshold. 

The exact effective continuum threshold — corresponding to exact values of the hadron mass and the decay constant 
on the left-hand side — is, of course, not known. Therefore, the actual extraction of hadron parameters from a sum rule 
consists in attempting (i) to find some reasonable approximation to the exact threshold and (ii) to control the accuracy 
of such an approximation. We stress again that the use of a r-dependent threshold is expected to improve the reliability 
of the extraction of the hadron parameter considered compared with the standard procedure of assuming a constant, 
r-independent threshold. 

Let us now look in detail at each step of the sum-rule calculation of the decay constant, starting with the OPE for 
the correlator. 



ds e 



Ppcrt ( 



oft(T) 



(1.7) 



2. OPE AND HEAVY-QUARK MASS 



We use the perturbative spectral density Pport(s) calculated in [3] to three-loop accuracy in terms of the pole 
mass of the heavy quark. The pole mass has been used in most of the sum-rule analyses since the pioneering work 
. An alternative option is to reorganize the perturbative expansion in terms of the running MS mass . Since the 
correlator is known to a^-accuracy, the relationship between pole and MS mass to the same accuracy is used. Explicit 
expressions for the perturbative spectral densities and power corrections may be found in 0, Q and are not given here. 

Figure [1] shows the perturbative spectral densities and the sum-rule estimates for Jd arising from (|1.8p for our two 
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choices of TOcI the pole mass mc,poic and the running MS mass 'mc{^)- The relevant OPE parameters are 

mc(mc) = (1.279 ± 0.013) GeV, m(2 GeV) = (3.5 ± 0.5) MeV, to,(2 GeV) = (100 ± 10) MeV, 

05(^2) = 0.1176 ± 0.0020, (2.1) 

{qq){2 GeV) = -((267 ± 17) MeV)^ (ss)(2 GeV)/(gg)(2 GeV) = 0.8 ± 0.3, ( — Gg) = (0.024 ± 0.012) GeV^ 



We employ a recent determination Q of mc{mc)- The corresponding pole mass recalculated from the 0{al) relation 
between Trie and TOc.poio is 



"ic,poic — 1.682 GeV. 



(2.2) 



The sum-rule estimates shown in Fig.[T]are obtained for a r-independent effective threshold Sq. Its values, which prove 
to be different for the pole-mass OPE and the MS-mass OPE, are found by requiring maximal stability of the extracted 
decay constant. Obviously, for heavy-light correlators and the resulting decay constants it makes a very big difference 
which precise scheme for the heavy-quark mass is employed. 



{ajn)' p,{s) 
0.4 

0.3 

0.2 

0.1 



-0.1 





Pole-mas^^^^ 




i=0i=li=2 Pull Pert 



f dud(T,so=4.7 GeV) 
0.04 

0.03 

0.02 

0.01 



Pole-mass 



0(1) 0(0-) 0(a^) power total 



(a/n)' p,(s,li=mc) 
0.4 



0.3 
0.2 
0.1 


-0.1 

al(T,So 

0.04 
0.03 
0.02 
0.01 









i=Oi=l i=2Full Pert 


2 3 4 5 e 

=5.0 GeV) 



MS-mass 



0(1) 0(ff) 0(a^) power total 



0.1 



0.2 0.3 0.4 0.5 



0.1 



0.2 



0.3 0.4 



0.5 



Fig. 1: OPE calculated in terms of the pole mass (left) and the MS mass (right) of the c quark. First line: spectral densities; 
second line: corresponding sum-rule estimates for fo- A constant effective continuum threshold so is fixed in each case separately 
by requiring "maximal stability" of the extracted decay constant. As the result, so turns out to be different in the two schemes. 



Several lessons should be learnt from these plots: 

(i) The perturbative expansion for the decay constant in terms of the pole mass shows no signs of convergence; each of 
the terms — LO, NLO, NNLO — gives contributions of similar size. Therefore, there is no reason to expect higher orders 
to give smaller contributions. 

(ii) Reorganizing the perturbative series in terms of the MS mass of the heavy quark leads to a clear hierarchy of the 
perturbative contributions. 

(iii) The absolute value of the decay constant extracted from the pole-mass OPE {fjj — 150 MeV) proves to be 
considerably smaller than that from the MS scheme {fjj — 180 MeV). Let us emphasize that, nevertheless, in both 
cases the decay constant exhibits perfect stability in a wide range of the Borel parameter r! Thus we emphasize again 
that mere Borel stability is by far not sufficient to guarantee the reliability of the sum-rule extraction of bound-state 
parameters. We have already observed this feature in several examples in quantum mechanics [6]. 

Because of the obvious problems with the pole-mass OPE for the correlator, we shall make use of the OPE in terms 
of the running MS mass for our extraction of the decay constants. Hereafter, the quark masses niQ and m, and the 
strong coupling as denote the MS running quantities at the scale ^. 
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3. EXTRACTION OF THE DECAY CONSTANT 

In order to determine the heavy-meson decay constant /g from the OPE, we must execute the following two steps. 

1. The Borel window 

First, we must fix our working r-window where, on the one hand, the OPE gives a sufficiently accurate description of 
the exact correlator (i.e., all higher-order radiative and power corrections are small) and, on the other hand, the ground 
state gives a "sizable" contribution to the correlator. Since the radiative corrections to the condensates increase rather 
fast with r, it is preferable to stay at the lowest possible values of r. We shall therefore fix the window by the following 
criteria (a) In the window, power corrections npowcr(T) should not exceed 30% of the dual correlator Ilduai{T, sq). 

This restricts the upper boundary of the r-window. The ground-state contribution to the correlator at this value of t 
comprises about 50% of the correlator, (b) The lower boundary of the r-window is fixed by the requirement that the 
ground-state contribution docs not fall below 10%. 

8. The effective continuum threshold 

Second, we must define a criterion how to determine Scs{t). The corresponding algorithm has been formulated in our 
recent works 0, and was shown to provide a good extraction of the ground-state parameters in quantum- mechanical 
potential models. 

Let us introduce the dual invariant mass iV/duai and the dual decay constant /dual by the relations 

MlAr) ^ lognduai(r, scff (r)), fl,,{T) ee MQ^e'^i^U^Ur, Seff (r)). (3.1) 

For a properly constructed Ilduai (''■, Scff (t) ) , this dual mass should coincide with the actual mass of the ground state. So, 
if the ground-state mass is known, any deviation of the dual mass from the actual mass of the ground state yields an 
indication of the contamination of the dual correlator by excited states. 

Assuming some particular functional form of the effective threshold and requiring the least deviation of the dual mass 
(13. ip from the actual mass in the r-window entails a variational solution for the effective threshold; as soon as the latter 
has been fixed, p.ip yields the decay constant. The standard assumption for the effective threshold is a r-independent 
constant. In addition to this approximation, we also consider polynomials in r. 

Our algorithm for the extraction of fg makes use of the knowledge of the true PQ-meson mass Mq. This algorithm, 
developed in our previous works and proven to work well for different correlators in the potential model, is very simple: 
we consider the set of r-dependent Ansatze for the effective continuum threshold 

n 

4ff^W=E4"^^^- (3-2) 

3=0 

We fix the parameters on the right-hand side of (13.21) as follows: we compute the dual mass squared according to p.ip 
for the r-dependent Scs{t) in (13. 2|) . We then evaluate M^^^^{t) at several values of r = r^ (i = 1, 2, . . . , A^, where A^ can 
be taken arbitrary large) chosen uniformly in the window. Finally, we minimize the squared difference between M^^^^ 
and the known value Mg: 

^'^]^E[^-^duai(T-0-M^]'. (3.3) 

1=1 

This gives us the coefficients s^"'' of the effective continuum threshold. As soon as the latter is fixed, it is straightforward 
to calculate the decay constant. 

The results presented below indicate that accounting for the r-dependence of the effective threshold yields a visible 
improvement compared with the usual assumption of a r-independent quantity in the following respect: it leads to a 
much better stability of the dual mass calculated for a dual correlator, which is tantamount to a better isolation of the 
ground-state contribution. 

Still, by trying different Ansatze for the effective continuum threshold, one obtains different estimates for the 
decay constant. We discuss the interpretation of these results in connection with the systematic uncertainties of the 
method of sum rules. 
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3. Uncertainties in the extracted decay constant 

Clearly, the extracted value of the decay constant is sensitive to the precise values of the OPE parameters and to the 
prescription for fixing the effective continuum threshold. The corresponding errors in the resulting decay constants are 
called the OPE-related error and the systematic error, respectively. Let us discuss these in turn. 

OPE-related error 

The value of the OPE-related error is obtained as follows: We perform a bootstrap analysis 9] by allowing the OPE 
parameters to vary over the ranges indicated in (|2.ip . using 1000 bootstrap events. Gaussian distributions for all OPE 
parameters but /i are employed. For /i we assume a uniform distribution in the corresponding range, which we choose 
to be 1 < /i (GeV) < 3 for charmed mesons and 2 < fi (GeV) < 8 for beauty mesons. The resulting distribution of the 
decay constant turns out to be close to Gaussian shape. Therefore, the quoted OPE-related error is a Gaussian error. 

Systematic error 

The systematic error of some hadron parameter determined by the method of sum rules (i.e., the error related to the 
intrinsic limited accuracy of this method) represents the perhaps most subtle point in the applications of this method. 
So far no way to arrive at a rigorous — in the mathematical sense — systematic error has been proposed. Therefore, in 
this respect we have to rely on our experience obtained from the examples where the exact hadron parameters may be 
calculated independently from the method of dispersive sum rules and then compared with the results of the sum-rule 
approach. Recent experience from potential models shows that the band of values obtained from linear, quadratic, and 
cubic Ansatze for the effective threshold encompasses the true value of the decay constant f?\. Moreover, we could show 
that the extraction procedures in quantum mechanics and in QCD are even quantitatively rather similar [8j] . Therefore, 
we believe that the half- width of this band may be regarded as realistic estimate for the systematic uncertainty of the 
decay constant. Presently, we do not see other possibilities to obtain a more reliable estimate for the systematic error. 



A. Decay constant of the D meson 



The r- window for the charmed mesons, r = (0.1-0.5) GeV-^ is chosen according to the criteria formulated above. 
Figure [2] shows the application of our procedure of fixing the effective continuum threshold and extracting the resulting 
fi). We would like to point out that, in the window, the r-dependent effective thresholds reproduce the meson mass 
much better than the constant one (Fig. [2^). This signals that the dual correlators corresponding to the r-dependent 
thresholds are less contaminated by excited states. 



Md„.^/MD fd„ai[MeV] fDlMeV] 




(a) (b) (C) 

Fig. 2: Dual mass (a) and dual decay constant (b) of the D meson obtained using different Ansatze for the effective continuum 
threshold Sch(t) 1)3. 2[) and fixing all thresholds according to (13. 3|) . Results for nic = mc{jnc) = 1.279 GeV, pi = mc, and central 
values of the other relevant parameters are presented, (c) Dual decay constant of the D meson vs. mc for /i — mc and central 
values of the other OPE parameters. The integer n = 0, 1, 2, 3 is the degree of the polynomial in our Ansatz (|3.2p for Scs{t). 



The dependence of the extracted value of the D-meson decay constant fu on the c-quark mass rUc 
condensate (qq) = {qq(2 GeV)) may be parameterized as 



mc{mc) and the 



j'dual 

Jd 



206.2 - 13 



1.279 GeV' 



0.1 GeV J 



I (gg) 1 - 0.267 GeV 
0.01 GeV 



±5.1 



(,y,t) MeV. (3.4) 



This formula describes the band of values indicated by the two dotted lines in Fig. which delimit the results found 
from the linear, quadratic, and cubic Ansatze for the effective continuum threshold. Figure[3ti. depicts the result of the 
bootstrap analysis of the OPE uncertainties. The distribution has a Gaussian shape, and therefore the corresponding 
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Fig. 3: (a) Distribution of fo obtained by the bootstrap analysis of the OPE uncertainties. Gaussian distributions for all OPE 
parameters but /x with corresponding errors as given in H2.1|l are employed. For u we assume a uniform distribution in the range 
1 GeV < n < 3 GeV. (b) Summary of findings for fo- Lattice results are from [13, EH for two dynamical light flavors {Nf = 2) 
and from [l^. [l3|] for three dynamical flavors {Nf = 3). The triangle represents the experimental value from PDG For the 
r-dependent QCD-SR result the error shown is the sum of the OPE and systematic uncertainties in (|3.5p . added in quadrature. 



OPE uncertainty is the Gaussian error. Adding the half- width of the band deduced from our r-dependent Ansatze for 
the effective continuum threshold of degree n — 1, 2, 3 as the (intrinsic) systematic error, we obtain the following result: 

fD = (206.2 ± 7.3(oPE) ± 5.1(syst)) MeV. (3.5) 

The main sources of the OPE uncertainty in the extracted fo are its renormalization-scale dependence and the error of 
the quark condensate. 

For a T-independent Ansatz for the effective continuum threshold a bootstrap analysis entails the substantially lower 
range 

j,(n=o) ^ (181.3 ±7.4 

(OPE)) MeV, which differs from our r-dependent result p.Sp by ~ 10%, i.e., by almost three 
times the OPE uncertainty. Moreover, as we have already shown in our previous works Q, making use of merely the 
constant Ansatz for the effective continuum threshold does not allow one to probe at all the intrinsic systematic error 
of the QCD sum rule. From our result (|3.5|) the latter turns out to be of the same order as the OPE uncertainty. 

Allowing the threshold to depend on t leads to a clearly visible effect and brings the results from QCD sum rules into 
perfect agreement with recent lattice results and the experimental data (Fig.[3lD). This perfect agreement of our result 
with both experimental data and lattice results provides a strong argument in favour of the reliability of our procedure. 

B. Decay constant of Da meson 

The corresponding r-window is t = (0.1 — 0.6) GeV^^. FigureUprovides the details of our extraction procedure. Our 



Mdud/Mo, fj„al[MeV] fDK[MeV] 




Fig. 4: Same as Fig. [2] but for the Ds meson. 



7 



results for the Ds-meson decay constant fo^ niay be represented as 



245.3- 18 



TUr - 1.279 GeV 



0.1 GcV / 



3.5 



1(55)1^^- 0.248 GeV \ 
0.01 GeV J 



±4.5 



(syst) 



MeV. (3.1 



This formula describes the band of values indicated by the two dotted lines in Fig. St as function of rric = mc{mc) and 
gives also the dependence on the quark condensate (ss) = (ss(2 GeV)). Performing the bootstrap analysis of the OPE 
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Fig. 5: (a) Distribution of fo^ obtained by the bootstrap analysis of the OPE uncertainties. Gaussian distributions for all OPE 
parameters but fi with corresponding errors as given in 1)2. ip are employed. For ^ we assume a uniform distribution in the range 
1 GeV < /i < 3 GeV. (b) Summary of findings for fo^ ■ Lattice results are from [13, [Ol for two dynamical light fiavors {Nf = 2) 
and from [13, [131 for three dynamical fiavors {Nj = 3). The triangle represents the experimental value from PDG 1^. For the 
r-dependent QCD-SR result the error shown is the sum of the OPE and systematic uncertainties in (|3.7p . added in quadrature. 



uncertainties, we obtain the following estimate; 

/d, = (245.3 ± 15.7(oPE) ± 4.5(,y,t)) MeV. 



(3.7) 



As in the case of /c a constant-threshold Ansatz yields a substantially lower value: 



(218.8 ± 16.1(oPE)) MeV. 



C. UJfD 

For the ratio of the D and Ds decay constants we report the sum-rule prediction 

Id J Id = 1.193 ± 0.025(ope) ± 0.007(,y,t). (3.8) 

This value is to be compared with the PDG average foj fo — 1.25 ±0.06 [l3| as well as with the recent lattice results 
foJfD = 1.24±0.03 [in for Nf = 2 and JdJ/d = l.i64±0.011 [H and fajfo = 1.20±0.02 [l| for TV/ 3. The 
error in (|3.8p arises mainly from the uncertainties in the quark condensates {ss)/{qq) = 0.8 ± 0.3. 



4. SUMMARY AND CONCLUSIONS 



We presented a detailed analysis of the decay constants of charmed heavy mesons with the help of QCD sum rules. 
Particular emphasis was laid on the study of the uncertainties in the extracted values of the decay constants: the OPE 
uncertainty related to the not precisely known QCD parameters and the intrinsic uncertainty of the sum-rule method 
related to a limited accuracy of the extraction procedure. 

Our main findings may be summarized as follows. 

(i) The perturbative expansion of the two-point function in terms of the pole mass of the heavy quark exhibits 
no sign of convergence. However, reorganizing this expansion in terms of the corresponding running mass leads to a 
clear hierarchy of the perturbative contributions. Interestingly, the decay constant extracted from the pole-mass OPE 
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proves to be sizeably smaller than the one extracted from the running-mass OPE. In spite of this numerical difference, 
the decay constants extracted from these two correlators exhibit perfect stability in the Borel parameter. This example 
shows that stability per se does not guarantee the reliability of the sum-rule extraction of any bound-state parameter. 

(ii) We have made use of the Borel-parameter-dependent effective threshold for the extraction of the decay constants. 
The r-dependence of the effective threshold emerges quite naturally when one attempts to increase the accuracy of the 
duality approximation. According to our algorithm, one should consider different polynomial Ansatze for the effective 
threshold and fix the coefficients in these Ansatze by minimizing the deviation of the dual mass from the known actual 
meson mass in the window. Then, the band of values corresponding to the linear, quadratic, and cubic Ansatze refiects 
the intrinsic uncertainty of the method of sum rules. The efficiency of this criterion has been tested before for several 
examples of quantum-mechanical models. This strategy has now been applied to the decay constants of heavy mesons. 

(iii) We obtained the following sum-rule estimates for the decay constants of the charmed D and Dg mesons: 

fo - (206.2 ±7.3(oPE)±5.1(,yst))MeV, (4.1) 
fo, = (245.3 ±15.7(oPE)±4.5(,y,t))MeV. (4.2) 

We point out that we provide both the OPE uncertainties and the intrinsic (systematic) uncertainty of the method of 
sum rules related to the limited accuracy of the extraction procedure. In the case of fn the latter turns out to be of the 
same order as the OPE uncertainty. Noteworthy, adopting a r-independent effective threshold leads to a substantially 
lower range = (181.3 ±7.4 

(OPE)) MeV, which differs from our r-dependent result (|4.1|) by almost three times 
the OPE uncertainty. The resulting ratio of the decay constants is 

IdJId - 1.193 ± 0.025(OPE) ± 0.007(,y,t)- (4.3) 

(iv) Our study of charmed mesons clearly demonstrates that the use of Borel-parameter-dependent thresholds leads 
to two essential improvements: 

a. The actual accuracy of the decay constants extracted from sum rules improves considerably. 

b. Our algorithm yields realistic (although not entirely rigorous) estimates for the systematic errors and allows one to 
reduce their values to the level of a few percent. Due to the application of our prescription, the QCD sum-rule results 
are brought into perfect agreement both with the experimental results and with lattice QCD. 
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